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Subject:  An  Experimental  Study  of  Travelling-Bubble  Cavitation  Noise 
References:  See  page  27. 

Abstract:  In  this  investigation,  high-frequency  travelling-bubble  cavitation 
noise  generated  on  a  Schiebe  headform  is  analyzed  with  regard 
to  noise  theories  which  are  based  on  the  collapse  of  a  single 
bubble.  After  normalizing  the  noise  level  by  an  event  rate, 

-o/C_  appears  to  provide  a  reasonable  scaling  factor  for 
*min 

the  noise  level  radiated  per  collapse,  where  o  is  the  cavitation 
number  and  C_  is  the  minimum  pressure  coefficient  on  the 
Hmin 

headform.  At  frequencies  above  10  kHz  and  for  free-stream 
flow  velocities  between  9.1  and  12.2  mps  (30  and  40  fps), 
the  noise  level  per  collapse  varied  as  roughly  0-8. 4 ^  and 

was  constant  with  the  flow  velocity  for  constant  -o/CD 

Fmin 

The  energy  density  spectra  were  flat  from  10  to  100  kHz 
regardless  of  the  cavitation  number  or  free-stream  velocity. 
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NOMENCLATURE 

radius  of  Schiebe  body- 

speed  of  sound 

pressure  coefficient 

average  pressure  coefficient  where  P  < 

minimum  pressure  coefficient 

frequency 

noise  level 

noise  level  per  collapse 
shock  wave  parameter  • 

*  '  j 

number  of  collapses  per  second 

acoustic  pressure 

pressure 

noncondensable  gas  pressure 

peak  shock  wave  pressure 

vapor  pressure 

free-stream  static  pressure 

radial  coordinate 

bubble  radius 

bubble  wall  velocity 

bubble  wall  acceleration 

minimum  bubble  radius 

maximum  bubble  radius 

length  on  surface  of  Schiebe  body 

length  of  region  where  P  < 

spectral  energy  density 

time 
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NOMENCLATURE  (Cont'd) 

bubble  volume 

second  derivative  of  bubble  volume 
free-stream  velocity 
time  constant 
mass  density  of  water 
cavitation  number 

9 

desinent  cavitation  number 
incipient  cavitation  number 
bubble  collapse  time  " 3 
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INTRODUCTION 

Although  there  have  been  many  experimental  investigations  of  cavitation 
noise,  most  studies  are  generally  concerned  with  well-developed  forms 
such  as  wake,  vortex  or  fixed-patch  cavitation.  In  those  cases,  the 
entire  cavitating  region  is  regarded  as  a  single  complex  noise  source 
as  opposed  to  considering  the  collapse  of  each  individual  bubble.  Near 
inception  when  the  bubbles  are  collapsing  independently,  the  radiated 
noise  is  often  more  intense  than  when  the  cavitation  becomes  developed. 

Such  incipient  cavitation  noise  has  already  been  analyzed  -  for  example 
on  hydrofoils  by  Barker  [1]  and  Blake  et.  al.  [2]  -  but  because  of  the 
typically  narrow  ranges  of  flow  parameters  which  support  incipient  conditions 
relations  indicating  how  the  noise  depends  on  the  flow  parameters  have 
not  been  determined  for  these  cases.  Noise  studies  have  also  been  conducted 
for  the  collapse  of  a  single  bubble.  This  was  accomplished  by  Harrison  [3], 
Mellen  [4],  and  more  recently  by  Chahine  et.  al.  [5]  where  spark  gaps  were 
used  to  generate  the  bubbles,  and  also  by  Hentschel  [6]  who  employed  a 
focused  laser  beam  for  that  purpose.  However,  the  parameters  involved  in 
those  Investigations  do  not  reflect  those  typically  encountered  in  flow- 
induced  cavitation. 

Therefore,  the  present  study  examines  travelling-bubble  cavitation 
which  is  easily  generated  on  a  Schiebe  headform.  In  travelling-bubble 
cavitation  each  bubble  collapses  independently,  while  the  Schiebe  headform 
allows  for  a  relatively  wide  operating  range,  in  terms  of  flow  parameters, 
for  this  type  of  cavitation  to  exist.  Measurements  of  cavitation  noise 
generated  in  this  way  are  compared  to  the  relevant  theories  for  noise 
radiation  due  to  bubble  collapse. 
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Cavitation  noise  analysis  essentially  began  in  1917  with  Rayleigh's 

solution  [7]  to  the  problem  posed  by  Besant  [8]  of  calculating  the  pressure 

field  generated  by  the  collapse  of  a  spherical  cavity  in  an  incompressible 

fluid.  The  extension  of  this  analysis  by  Plesset  [9]  resulted  in  the 

Rayleigh-Plesset  equation  governing  the  dynamics  of  a  bubble  in  an 

arbitrary  time-varying  pressure  field.  It  has  been  shown  [10]  how  this 

classical  analysis  predicts  an  energy  spectrum  which  rises  as  roughly 

4 

f  to  a  peak  at  a  frequency  near  the  reciprocal  of  the  bubble  lifetime, 

-2/5 

after  which  it  rolls  off  as  f  .  Both  Mellen  [11]  and  Blake  et.  al.  [2] 

obtained  noise  spectra  which  exhibited  this  type  of  spectral  response, 

-2 

except  that  their  roll-offs  behaved  as  f 

Subsequent  theoretical  advances  in  modelling  cavitation  noise  were 
few  and  far  between.  Not  until  World  War  II  was  an  analysis  of  the  propa¬ 
gation  of  weak  shocks  developed  by  Kirkwood  and  Bethe  [12]  which  enabled 
Gilmore  [13]  in  1952  to  derive  an  equation  for  the  collapse  of  a  bubble 
in  a  compressible  fluid.  Both  experiment  [4]  and  numerical  analysis 
[14]  have  demonstrated  the  accuracy  of  Gilmore's  theory,  particularly 
around  the  final  stage  of  collapse  where  bubble  wall  velocities  are  sonic 
and  the  effects  of  compressibility  become  important.  However,  solutions 
of  that  equation  must  be  obtained  by  numerical  methods,  and  as  such  it 
does  not  yield  useful  analytical  predictions  of  the  radiated  sound  field. 

Another  important  parameter  during  the  final  collapse  stage  is  the 
presence  of  any  noncondensable  gas  contained  in  the  bubble.  In  addition, 
acoustic  pressure  waves  may  be  generated  by  the  high  bubble  wall  velocities 
Using  Whitham's  method  [15]  of  describing  shock  wave  propagation.  Baiter 
[16]  recently  derived  an  expression  for  the  pressure  pulse  in  the  far 
field  of  the  bubble,  including  the  effects  of  noncondensable  gas  pressure 
but  not  of  compressibility  with  respect  to  bubble  dynamics.  Esipov  and 
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Naugol'nykh  [17],  using  the  Kirkwood-Bethe  approximation,  did  include 
the  effects  of  compressibility  and  arrived  at  a  solution  similar  to  Baiter's 
Both  solutions  assume  a  shock  wave  having  an  exponential  profile,  and 
consequently  depict  a  flat  energy  spectrum  which  eventually  rolls  off 
at  the  high  end  around  a  frequency  determined  by  the  time  constant 
of  the  wave  profile.  A  number  of  experiments  [1,2,11]  yielded  noise 
spectra  which  hinted  at  such  a  high  frequency  spectral  plateau. 

Thus,  this  investigation  will  relate  these  noise  theories  to  the 
sound  radiated  by  travelling-bubble  cavitation.  Comparisons  are  based 
on  how  the  noise  level  per  collapsing  bubble  depends  on  the  cavitation 
number  a  and  the  flow  velocity  V  .  The  independent  behavior  of  the  bubbles 
permits  the  measured  frequency  spectra  to  be  compared  with  those  predicted 
by  the  single  bubble  noise  theories  [18]. 

EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

This  investigation  was  conducted  in  the  305  mm  (12  in.)  diameter 
water  tunnel  of  the  Applied  Research  Laboratory  at  the  Pennsylvania  State 
University.  The  headforms  used  were  two  50.8  mm  (2  in.)  diameter  Schiebe 
noses  which  were  identical  except  for  modifications  that  provided  measure¬ 
ments  of  either  static  or  acoustic  pressure.  Schiebe  headforms  do  not 
experience  laminar  separation,  and  they  support  travelling-bubble  cavitation 
for  a  relatively  wide  pressure  range  provided  there  are  enough  microbubbles 
available  in  the  free  stream  [19]. 

To  measure  the  pressure  distribution  along  the  surface  of  the  Schiebe 

nose,  one  headform  was  equipped  with  ten  flush-mounted  pressure  taps 

at  various  distances  from  the  stagnation  point.  These  pressures  were 

compared  with  the  free-stream  static  pressure,  for  flow  velocities  ranging 

from  6.1  to  18.3  raps  (20  to  60  fps),  to  determine  the  pressure  coefficient 

C  defined  as 
P 


gggai 


-9- 


4  December  1981 
MFH : DET : MLB : cag 


With  the  cavitation  number  a  given  by 

P  -  P 


oo  v 


(l) 


(2) 


measurements  of  the  incipient  and  desinent  cavitation  numbers,  and 
od  respectively,  were  determined  visually  with  the  aid  of  a  stroboscope 
by  adjusting  the  test  section  pressure  at  constant  flow  velocities. 

Attempts  to  use  the  cavitation  noise  to  define  these  quantities  are  described 
in  Reference  [20],  but  those  approaches  were  unreliable. 

A  study  of  bubble  dynamics  was  accomplished  with  a  high-speed  photographic 
system  which  consisted  of  a  Redlake  Hycam  camera  and  an  EG  &  G  type  501 
stroboscope.  The  Hycam  is  a  16  mm  rotating  prism  camera  which  was  adjusted 
to  a  framing  rate  of  roughly  1800  pictures  per  second.  A  Schiebe  nose 
was  modified  with  a  small  barium- ti tana te  transducer  which  was  mounted 
flush  with  the  surface  of  the  model  60  mm  downstream  from  the  stagnation 
point.  Thus,  at  the  same  time  that  the  high-speed  films  of  bubble  collapse 
were  taken,  the  resulting  noise  was  being  monitored  by  that  transducer. 

After  high-pass  filtering  at  10  kHz,  the  noise  was  recorded  on  a  Bell  and 
Howell  type  CPR-4010  tape  recorder.  The  tape  recording  was  synchronized 
with  the  films  via  voltage  spikes  input  to  a  second  channel  which  marked 
the  points  in  time  where  the  strobe  would  fire  to  create  each  successive 
picture. 

Additional  studies  of  bubble  dynamics  were  realized  by  filming  the 
cavitation  on  videotape.  Voltages  corresponding  to  the  pressure  and 
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velocity  In  the  test  section  were  called  out  and  recorded  on  the  audio 
track  of  the  tape  so  that  the  effects  of  a  wide  range  of  flow  parameters 
could  be  observed. 

To  measure  the  cavitation  noise  level,  a  25  mm  (1  In.)  diameter 
lead-zirconate-titanate  (PZT)  hydrophone  was  coupled  to  the  outside  of 
the  test  section  window.  The  specific  acoustical  impedance  of  the  window 
was  near  that  of  water  such  that  minimum  sound  reflection  occurred  at 
the  interface  formed  by  the  two  media.  The  hydrophone  was  aimed  at  a 
location  on  the  model  about  50  mm  downstream  from  the  stagnation  point. 

With  a  separation  of  about  150  mm  (6  in.)  between  the  face  of  the  hydrophone 

•  *  5 

and  the  surface  of  the  model,  the  measured  beamwidth  on  that  surface 
was  about  70  mm  at  its  half-power  points. 

After  high-pass  filtering  at  10  kHz,  the  output  from  the  hydrophone 
was  input  to  one  channel  of  the  tape  recorder.  A  DC  voltage  output  from 
an  rms  voltmeter,  which  represented  the  cavitation  noise  level  above 
10  kHz,  was  input  to  a  second  channel.  On  yet  a  third  channel  was  recorded 
a  DC  voltage  which  was  proportional  to  the  free  stream  static  pressure 
in  the  test  section.  The  entire  arrangement  is  illustrated  in  Figure  1. 

Tests  were  run  at  flow  velocities  ranging  from  9.1  to  12.2  mps 
(30  to  40  fps)  in  increments  of  0.76  mps  (2.5  fps) .  With  the  velocity 
held  constant,  the  pressure  in  the  test  section  was  lowered  from  an 
initial  point  above  o^..  As  soon  as  noise  bursts  due  to  bubble  collapse 
were  witnessed  on  an  oscilloscope,  tape  recording  would  begin  as  the 
pressure  was  continuously  lowered  until  the  travelling-bubble  cavitation 
transformed  into  an  attached  cavity.  This  procedure  was  repeated  at  each 
velocity.  In  addition,  noise  spectra  from  10  to  100  kHz  were  obtained 


-11- 


4  December  1981 
MFH : DET : MLB : cag 


using  a  Celesco  LC-10  hydrophone  mounted  In  a  parabolic  reflection  dish 
on  the  test  section  window.  Adjustments  for  natural  resonances  due  to 
the  test  section  geometry  were  made  on  the  basis  of  spectra  generated  by 
a  calibrated  LC-10  source  located  adjacent  to  the  model.  The  spectra 
were  analyzed  with  a  Spectral  Dynamics  SD-360  real  time  FFT  processor. 

Finally,  the  effect  of  the  velocity  and  pressure  on  the  sound 
transmission  path  was  investigated.  The  temperature  of  the  water  used 
throughout  this  investigation  was  approximately  21°C,  while  the  total 
gas  content  was  roughly  10  parts  per  million  (PPM)  air  to  water  on  a 
molar  basis  as  measured  with  a  van  Slyke  apparatus.  To  determine  how  sound 
was  absorbed  by  free-stream  bubbles  as  a  function  of  pressure,  white 
noise  was  projected  into  the  test  section  through  the  window  opposite 
the  PZT  hydrophone.  This  test  was  conducted  at  7.6  mps  (25  fps)  since 
lower  pressures  are  required  to  support  cavitation  at  correspondingly  lower 
velocities,  and  hence  any  effects  due  to  free-stream  bubbles  should  be 
prevalent  under  these  conditions.  Attenuation  of  less  than  1  dB  was 
observed  above  10  kHz  after  lowering  the  pressure  from  a  state  of  no 
cavitation  until  the  cavitation  noise  itself  obscured  the  projected  white 
noise.  Consequently,  this  effect  was  ignored. 

There  were  essentially  two  reasons  for  high-pass  filtering  all  noise 
measurements  at  10  kHz.  First,  the  water  tunnel  itself  generated  so  much 
noise  at  low  frequencies  as  to  render  making  adjustments  for  the  background 
noise  almost  impossible.  Second,  near-field  effects  make  interpretation 
of  measurements  below  10  kHz  questionable.  The  wavelength  of  sound  in 
water  at  that  frequency  is  already  about  150  mm  (6  in.),  which  is  approxi¬ 
mately  the  distance  from  the  region  of  cavitation  to  the  hydrophone. 
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MATHEMATICAL  MODELS 

The  dynamics  of  a  spherical  bubble  in  an  incompressible  fluid  are 
described  by  the  Rayleigh-Plesset  equation  [9], 


••  3*2 

R  R  +  |  R 


P(R)  -  P_ 


where  P(R)  represents  the  pressure  applied  by  the  contents  of  the  bubble 
to  the  bubble  wall  located  at  radius  R  »  R(t),  and  P^  is  the  pressure 
field  experienced  by  the  bubble.  Very  often  P(R)  is  considered  to  be 
constant  and  equal  to  the  vapor  pressure  of  water,  P  .  This  approximation 
becomes  unacceptable  only  at  the  very  beginning  and  end  of  the  bubble 
life  cycle  where  gas  content  must  be  considered  for  small  bubble  volumes. 
Also,  effects  of  viscosity  and  surface  tension  have  been  shown  to  exert 
only  a  minor  influence  on  bubble  dynamics  [21]. 

Equation  (3)  yields  several  analytical  expressions  describing  various 
stages  of  growth  and  collapse  for  a  bubble  In  a  local  pressure  field 
defined  by  P.  During  vaporous  bubble  growth.  It  can  be  shown  [22]  that 
the  bubble  wall  velocity  remains  essentially  constant.  Thus,  setting 
R  “  0  and  P(R)  *  Py  in  Equation  (3),  an  approximate  formula  for  the  maximum 
bubble  radius  R^  is  easily  derived  and  given  by 


where  s ;  is  the  length  of  the  bubble  path  along  which  P  <  P  .  It  is 
assumed  here  that  bubble  growth  occurs  only  in  that  region,  with  the 
average  pressure  represented  by  C^. 

Even  during  collapse,  the  noncondensable  gas  pressure  becomes  significant 
only  for  R  <<  R^.  However,  at  this  point  the  bubble  wall  velocity  is  so 
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high  that  It  violates  the  assumption  of  Incompressibility,  invalidating 

Rayleigh's  analysis.  Nevertheless,  Rayleigh  calculated  the  collapse 

time  of  a  bubble  for  P(R)  =  P  to  be 

v 


(5) 


Hickling  and  Plesset  [14]  have  shown  Equation  (5)  to  agree  within  1Z  of 
their  numerical  analysis  for  a  compressible  liquid  in  the  absence  of  any 
vapor  pressure. 

In  the  far  field  of  the  bubble,  the  Rayleigh-Plesset  equation  yields 
input  for  calculating  the  acoustic  pressure  [23] 


p(r>t)  "  4^V(t  "f>  *  (6) 

3 

where  V  =  (4/3)  ttR  is  the  bubble  volume.  Using  the  Rayleigh-Plesset  equation 
again  to  derive  an  asymptotic  expression  for  the  radius  R(t)  of  a  collapsing 
bubble  for  R  «  R^,  the  acoustic  pressure  becomes 


P(r,t)  a  £ 


v„6/5  "h975 


(c  +<03/5t'4/5 
P 


Taking  the  Fourier  transform  of  Equation  (7)  and  squaring  the  result,  the 
energy  density  is  then 


S(f)  a  * 


v.12/5  <C 


+  „)6'5  f-2/5 


Any  noncondensable  gas  contained  in  a  bubble  cushions  the  collapse 
and  for  small  amounts  of  gas  the  collapsing  bubble  can  rebound  into  another 
growth  stage.  This  requires  an  abrupt  change  in  the  movement  of  the 
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surrounding  fluid  and  produces  a  compression  wave.  This  compression  wave 
propagates  outwards  and  can  possibly  form  a  shock  front.  Shock  formation 
is  caused  by  nonlinear  compressibility  effects  on  sound  radiation.  It 
occurs  when  the  initial  pulse  amplitude  is  high  enough  so  that  the 
shock  is  formed  before  the  nonlinear  effects  are  neutralized  through 
geometric  attenuation  and  propagation  losses.  Assuming  that  an  exponentially 
shaped  shock  wave  occurs.  Baiter  [16]  followed  Whitham's  analysis  [15]  and 
derived 

p(r,t)  -  P  e"t/0  ,  (9) 

s  ••  j 


to  arrive,  via  the  Rayleigh-Plesset  equation,  at  values  of  P  and  0  given 

8 

by 


P 

8 


a  pc 


2  *M 


(10) 


and 


(ID 


R  is  the  minimum  radius  of  the  bubble  at  the  instant  where  rebound  occurs, 
m 

and  Pg  refers  to  the  gas  pressure  in  the  bubble  when  R  -  R^j.  The  above 
relations  require  P  <<  P  for  a  gas  behaving  adiabatically  with  a  specific 
heat  ratio  equal  to  4/3. 

Considering  again  a  noncondensable  adiabatic  gas  but  accounting 
for  compressibility  of  the  fluid,  Esipov  and  Naugol'nykh  [17]  used  the 
Kirkwood-Bethe  approximation  [12]  to  arrive  at  an  expression  in  the  form 
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of  Equation  (9)  where  the  pressure  is  again  given  by  Equation  (10) ,  but 


—  ~ 


They  define  the  parameter  M  as 


M=  | 


3p 3/4  (pc2)i/4  r1 


Shock  wave  pulses  are  often  modelled  by  an  exponential  pulse.  The 

.  *  *  3 

energy  density  spectrum  of  Equation  (9)  is 


<v> 

1  4-  (2TV0f) 


2 


Equation  (14)  implies  a  flat  spectrum  up  until  a  cutoff  frequency  given  by 


after  which  the  spectrum  will  roll  off  as  f  .  Thus,  to  find  the  energy 
density  below  the  cutoff  frequency,  it  is  sufficient  to  determine 


S(f)  *»  (Pg0)  . 


For  the  range  of  flow  parameters  used  in  this  investigation,  the 

bubbles  experience  pressures  between  0.2  and  0.5  atm  during  their  collapse 

-3/2 

stage.  Consequently,  the  M  term  in  Equation  (12)  may  be  neglected 
for  P  <  0.001  atm.  In  fact,  such  low  gas  pressure  is  required  for  using 

O 
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Equations  (9)  through  (12)  because  of  the  assumption  that  shock  formation 
does  indeed  occur.  Hickling  and  Plesset  (14)  showed  that,  with  a  pressure 
of  1  atm  forcing  collapse,  shock  formation  just  barely  occurs  for  an 
initial  gas  pressure  of  0.01  atm.  Higher  gas  pressures  will  cushion 
the  collapse  sufficiently  to  prohibit  shock  formation  altogether.  Thus, 
Equation  (12)  becomes  almost  identical  to  Equation  (11)  so  that  the 
energy  density  below  cutoff  becomes,  from  Equation  (16), 

s<f)  «  *H  Pg  *  (17) 

r  6 

RESULTS  AND  DISCUSSION 

The  shape  of  a  Schiebe  headform  is  determined  by  the  profile  of 
streamlines  resulting  from  the  addition  of  a  disk  source  to  a  uniform 
flow  field  [24].  As  such,  no  closed-form  solution  describing  the  pressure 
distribution  along  its  surface  exists,  so  a  computer-generated  solution 
appears  in  Figure  2.  Distance  along  the  surface  of  the  headform  is 
measured  by  the  dimensionless  quantity  s/a,  where  the  actual  length  s 
is  normalized  by  the  maximum  radius  of  the  model.  In  this  case,  that 
radius  is  a  °  25.4  mm  (1  in.).  The  experimental  values  of  appearing 
in  Figure  2  result  from  averaging  a  number  of  runs  through  the  range  of 
flow  velocities.  An  uncertainty  of  ±0.01  is  associated  with  each  point. 

Figure  3  gives  values  obtained  for  the  incipient  and  desinent  cavitation 
numbers,  and  respectively,  with  an  uncertainty  of  ±0.03.  The  error 
is  often  larger  at  lower  velocities  where  the  free-stream  pressure  must 
correspondingly  be  lower  to  cause  cavitation.  At  those  lower  pressures, 
free-stream  air  bubbles  can  become  large  enough  to  be  confused  with  cavitation 


on  the  model. 
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The  point  at  which  travelling-bubble  cavitation  abruptly  transforms 

into  an  attached  cavity  is  denoted  in  Figure  3  by  Attached  Cavity  a 

Note  that  this  parameter  remains  relatively  constant  for  different  velocities 

Thus,  as  the  velocity  Increases,  the  operating  range  of  pressure  that 

can  support  travelling-bubble  cavitation  diminishes.  This  established 

12.2  mps  (40  fps)  as  a  practical  upper  limit  for  the  flow  velocity  where 

meaningful  trends  of  noise  versus  pressure  could  be  measured. 

An  interesting  feature  of  Figure  3  results  from  a  comparison  with 

Figure  2.  Classical  theory  predicts  a.  “  o ,  =  -C  for  vaporous  cavitation 

Pmin 

implying  that  vaporous  cavitation  should  not  exist  for  any  a  >  -C  as 

••  >  Pmin 

there  is  then  no  location  on  the  model  where  P  <  P  .  With  Cn  =  -1.0 

~  v  Pmin 

on  the  headform,  an  alternative  type  of  cavitation  above  a  =  1.0  must 

be  suggested  to  account  for  the  values  of  a.  and  a,  at  velocities  below 

1  d 

11.4  mps. 

The  most  probable  cause  for  this  discrepancy  is  that  the  cavitation 
nuclei  contain  gas  having  pressure  greater  than  Py.  Thus,  passage  through 
a  region  of  local  pressure  less  than  that  of  the  gas  pressure  may  cause 
the  nuclei  to  grow  enough  to  be  confused  with  vaporous  cavitation.  For 
example,  o  *  1.2  at  a  velocity  of  9.1  mps  Indicates  a  minimum  pressure  of 
about  0.1  atm  on  the  surface  of  the  headform,  implying  a  gas  pressure 
of  that  same  value  as  compared  with  P  =0.03  atm.  As  the  bubble  grows, 

_3 

however,  the  gas  pressure  decreases  as  R  for  isothermal  expansion  or  as 
-4 

R  if  the  expansion  is  adiabatic.  This  results  from  gas  diffusion  through 
the  bubble  wall  being  insignificant  in  relation  to  the  rate  of  bubble 
growth  and  collapse.  Consequently,  the  pressure  rapidly  approaches  Py 
and  remains  essentially  constant  during  bubble  growth.  Thus,  vaporous 
cavitation  will  dominate  as  the  pressure  is  lowered. 
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The  trend  for  o^  and  to  decrease  with  increasing  velocity  is 
not  uncommon  at  high  air  contents  such  as  10  ppm.  Since  this  phenomenon 
is  discussed  more  thoroughly  in  Reference  [25],  it  is  sufficient  here  to 
note  that  the  effect  of  any  gas  content  is  to  increase  the  point  of  inception 
for  decreasing  flow  velocities.  Where  a  specific  example  at  9.1  mps  was 
shown  to  suggest  a  gas  pressure  of  0.1  atm,  this  merely  set  an  upper  bound 
for  the  gas  pressure  in  those  particular  cavitation  nuclei  which  were 
prone  to  expansion  at  that  value  of  a.  Since  the  nuclei  possess  a  distri¬ 
bution  in  size,  those  participating  in  cavitation  inception  represent 
the  largest  nuclei  in  the  spectrum.  Smaller  free-stream  nuclei  may  have 
lower  gas  pressures  because  of  higher  surface  tension  which  resists  gas 
diffusion  through  the  bubble  wall. 

Data  resulting  from  the  high-speed  films  and  videotapes  were  compared 
with  the  equations  describing  bubble  growth  and  collapse.  Measurements 
of  exceeded  the  predictions  of  Equation  [4]  by  factors  ranging  from 
2  to  4.  It  is  interesting  that  Blake  et.  al.  [2]  performed  similar  calculations 
for  travelling-bubble  cavitation  on  a  hydrofoil  and  found  Equation  (4) 
to  overestimate  their  maximum  bubble  radii  by  about  20%.  At  the  very 
least,  this  suggests  that  Equation  (4)  is  more  suited  to  the  experiment 
conducted  by  Blake  et.  al.,  therefore  inviting  comparison. 

A  critical  parameter  in  Equation  (4)  is  the  distance  sy  travelled  by 
the  bubble  where  it  experiences  a  local  pressure  P  <  P^.  The  illustrative 
examples  discussed  in  Blake's  investigation  [2]  show  sy  to  be  greater 
than  30  mm,  such  that  the  bubbles  have  almost  attained  their  maximum 
radii  by  the  time  they  leave  that  low  pressure  region.  For  the  pressure 
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range  encountered  In  the  present  investigation,  this  low  pressure  region 
extends  for  less  than  4  mm  along  the  Schiebe  headform.  Also,  the  high¬ 
speed  films  showed  that  bubble  growth  did  not  usually  begin  until  about 
the  time  that  the  nucleus  was  leaving  the  low  pressure  region  and  most 
of  the  bubbles  occurred  outside  of  the  boundary  layer  on  the  Schiebe 
headform. 

Thus  it  is  not  surprising  that  the  bubble  growth  observed  by  Blake 
et.  al.  is  adequately  modelled  by  Equation  (4),  where  the  growth  stage 
was  assumed  to  be  entirely  limited  to  the  region  where  P  <  Pv>  Yet 
clearly,  inertial  effects  responding  to  thp  low  pressure  impulse  are 
primarily  responsible  for  growth  on  the  Schiebe  headform.  What  is  not 
certain  is  whether  Equation  (4)  correctly  describes  the  functional  dependence 
of  on  o  and  Cp  despite  an  inability  to  predict  specific  values. 

Similarly,  measurements  of  collapse  times  were  compared  with  the 
predictions  of  Equation  (5).  The  measured  times  were  generally  longer 
than  those  calculated  using  Equation  (5)  by  a  factor  of  2  to  3,  while 
Blake  et.  al.  observed  collapse  times  only  15%  in  excess  of  the  predicted 
values.  Air  content  is  one  basic  difference  here;  Blake  et.  al.  held  the 
air  content  below  2.5  ppm,  but  in  this  investigation  it  was  approximately 
10  ppm.  Higher  air  content  yields  an  increased  noncondensable  gas  content 
which  will  retard  the  rate  of  collapse,  while  Equation  (6)  assumes  the 
pressure  in  the  bubble  to  be  constant  and  equal  to  Pv>  It  is  interesting 
that  Harrison  [3],  observing  bubble  cavitation  in  a  venturi  nozzle  where 
he  maintained  an  air  content  of  about  5  ppm,  recorded  data  for  bubble 
collapse  corresponding  exactly  to  Rayleigh's  theory. 
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A  second  difference  is  the  shape  of  the  pressure  curve  on  the  body. 
The  Schiebe  headform  has  a  very  steep  rise  in  pressure  downstream  of  the 
location  of  the  minimum  pressure.  This  can  be  contrasted  to  a  gentle  rise 
in  pressure  downstream  of  the  location  of  the  minimum  pressure  on  the 
hydrofoil  used  by  Blake  et.  al.  Equation  (5)  estimates  the  collapse  time 
on  the  assumption  of  an  averaged  constant  pressure;  but  to  be  realistic, 
one  must  allow  for  a  changing  pressure. 

As  verified  repeatedly  in  previous  investigations,  the  high-speed 
films  showed  sound  radiation  occurring  at  the  final  instant  of  collapse. 
Only  a  few  rebounds  were  observed  at  7.6  'fops  (25  fps),  while  none  were 
directly  observed  at  higher  velocities.  Although  those  runs  were  all 
performed  at  an  air  content  of  10  ppm,  additional  films  were  obtained  at 
7.6  mps  with  an  air  content  of  5  ppm.  Under  these  conditions  rebounding 
was  frequent  with  significant  sound  generated  by  the  collapse  of  the 
rebound.  In  fact,  Chahine  et.  al.  [5]  observed  pressure  pulses  due  to 
the  collapse  of  rebounds  that  were  as  large  and  sometimes  even  larger 
than  the  pulse  from  the  initial  collapse.  The  results  they  cited 
were  for  spark-induced  bubbles  collapsing  at  a  distance  less  than  the 
maximum  bubble  diameter  away  from  a  rigid  wall.  Cavitation  on  headforms 
occurs  well  within  that  proximity,  so  that  cavitation  noise  theories 
should  account  for  this  phenomenon  once  rebounding  becomes  significant. 

The  noise  radiation  per  bubble  collapse  was  analyzed  by  digitizing 
tape  recordings  of  the  cavitation  noise  and  viewing  the  output  on  an 
oscilloscope.  Digitization  was  performed  at  an  effective  sampling  rate 
of  200  kHz  such  that  the  acoustic  signal  was  analyzed  for  continuous 
intervals  of  185  msec  in  real  time.  For  each  time  segment,  the  noise 
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level  and  pressure  in  the  test  section  were  determined  by  voltages  recorded 
on  separate  channels,  which  were  sampled  immediately  before  and  after 
every  185  msec  interval.  Thus,  by  counting  the  number  of  noise  bursts 
occurring  in  each  segment,  the  manner  in  which  the  noise  level  per  collapse 
depended  on  the  pressure  was  determined.  This  procedure  was  repeated 
at  different  velocities. 

It  should  be  noted  that  this  measured  event  rate  is  not  necessarily 
the  same  as  the  number  of  cavitation  bubbles  being  created  at  that  time. 
However,  acoustic  reflections  off  the  tunnel  walls  were  found  to  be  insig¬ 
nificant  by  observing  the  response  of  acoustic  impulses  generated  in 
the  test  section.  Also,  since  the  noise  tests  were  conducted  at  velocities 
between  9.1  and  12.2  mps  where  the  high-speed  films  and  videotapes  indicated 
negligible  rebounding  action,  this  factor  was  ignored.  Finally,  the 
collapse  of  some  bubbles  may  be  accompanied  by  an  insignificant  sound 
output,  but  these  cases  were  assumed  to  represent  anomalous  behavior. 

Thus  it  may  be  assumed  that  the  number  of  noise  bursts  adequately  represents 
the  actual  number  of  cavitation  bubbles.  More  importantly,  good  relations 
for  the  noise  per  collapse  are  obtained  which  can  then  be  compared  with 
the  relevant  theories.  In  what  follows,  it  is  assumed  that  each  noise 
burst  counted  represents  a  single  collapse. 

Accordingly,  the  number  of  collapses  per  second,  N,  appears 
in  Figure  4  as  a  function  of  pressure  and  velocity.  The  cavitation  number 
o  is  normalized  by  o^  and  plotted  on  a  logarithmic  scale.  Lines  drawn 
through  the  data  at  each  velocity  to  fit  the  majority  of  points  in  a  mean 
square  sense  yield  slopes  which  are  quite  similar,  having  as  an  average 
functional  relationship  N  a  Similarly,  Figure  5  displays 

the  data  representing  the  dependence  of  the  total  noise  level,  L,  on 
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pressure  and  velocity.  Again  the  lines  drawn  through  the  data  possess 

-19  0 

similar  slopes,  which  yield  the  average  relationship  Lao  .  The 
uncertainty  in  the  exponents  is  ±1.8. 

Thus,  although  the  noise  level  increases  rapidly  as  the  pressure 
is  lowered,  so  does  the  rate  of  bubble  collapse.  Figure  5  shows  that 
the  noise  tends  toward  a  saturation  level  for  very  low  pressures.  Such 
a  trend  is  also  observed  to  a  lesser  extent  in  Figure  4  where  the  collapse 
rate  appears  to  level  off  shortly  before  the  abrupt  transition  into  an 
attached  cavity. 

Although  normalization  by  a ^  in  Figures  4  and  5  spread  out  the  data 
for  clearer  presentation,  it  was  initially  thought  that  o/o,,  or  equivalently 

•  •  3  Q 

a/a^,  would  serve  as  an  appropriate  scaling  factor  for  the  noise  per 

collapse.  The  reason  that  this  is  not  the  case  relates  to  the  discussion 

of  why  and  are  not  constant  in  Figure  3  where  predominantly  gaseous 

cavitation  was  assumed  to  exist  above  a  *  1.0  at  all  velocities.  This 

does  not  mean  that  nuclei  with  gas  pressures  as  high  as  0.1  atm  will 

not  give  rise  to  bubbles  that  will  eventually  collapse  with  sufficient 

violence  to  generate  appreciable  sound.  After  all,  for  adiabatic  expansion 

the  nuclei  need  only  grow  to  a  maximum  radius  just  over  three  times  their 

initial  radius  to  attain  =  0.001  atm  at  R  =  R^,  such  that  shock  formation 

is  almost  inevitable.  However,  since  gaseous  bubble  growth  is  not  nearly 

as  explosive  as  that  of  vaporous  cavitation,  cavitation  above  a  =  1.0 

may  indeed  produce  negligible  sound  output.  In  view  of  this,  it  is  not 

surprising  that  the  noise  does  not  scale  with  a,.  Instead,  -C  was 

d  p  . 

mm 

chosen  to  be  the  normalization  factor  since  it  corresponds  in  theory 
to  the  values  of  c  marking  vaporous  cavitation  inception. 
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Therefore,  the  curves  in  Figures  4  and  5  are  plotted  again  in 

Figures  6  and  7  against  -o/C  ,  where  -C  seems  a  reasonable 

Pmin  Pmin 

choice  for  the  normalization  factor.  In  fact,  the  curves  for  the  total 

noise  level  in  Figure  7  appear  to  have  collapsed  quite  well.  The  somewhat 

greater  scatter  associated  with  the  curves  in  Figure  6  may  simply  be  a 

reflection  of  the  difficulties  inherent  in  counting  the  noise  bursts, 

especially  with  increasing  cavitation  activity. 

The  choice  of  -C  for  the  normalizing  factor  is  based  on  the 
Pmin 

theory  that  vaporous  cavitation  begins  at  values  of  a  corresponding  to 

-C  .  To  further  verify  this  choice  of  normalizing  factors,  a  series 
Pmin 

of  tests  should  be  conducted,  following  the  present  scenario,  with  a  series 

of  bodies  having  different  values  of  C  .  Correlation  of  the  results 

Pmin 

using  -C  would  then  establish  its  validity  as  the  proper  normalizing 
Pmin 

factor. 

To  derive  relationships  for  the  noise  level  per  collapse,  Lc,  the 

curves  fitting  the  data  for  each  velocity  in  Figure  7  were  normalized 

by  the  corresponding  curves  in  Figure  6.  The  resulting  curves  are  displayed 

-8  4 

in  Figure  8  where  the  slopes  suggest  the  average  relationship  Lc  a  a  ", 
with  an  uncertainty  of  ±2.5  in  the  exponent.  The  salient  feature  of 
Figure  8  is  that  the  noise  per  collapse  appears  to  scale  with  -o/C 

p  . 
min 

alone,  with  apparently  no  contribution  from  V  . '  Thus,  for  each  free- 
stream  velocity,  the  functional  dependence  of  the  noise  on  -o/C 

Pmin 

is  roughly  the  same.  Additionally,  the  dependence  of  o  itself  on  V 

CO 

is  such  that  for  a  constant  value  of  -o/C  ,  the  noise  per  collapse 

Pmin 

is  practically  equivalent  for  any  value  of  V  . 

Aside  from  difficulties  in  counting  the  noise  bursts,  another  possible 
reason  for  the  scatter  of  the  slopes  in  Figure  8  is  that  the  conditions 


in  the  tunnel  change  slightly  from  one  run  to  the  next.  This  was  borne 
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out  by  the  different  values  of  both  a^  and  when  measured  at  different 
times.  An  important  parameter  affecting  those  results  is  the  average 
size  of  the  free-stream  nuclei  which  is  known  to  change  significantly 
even  during  a  single  test.  Preliminary  data  obtained  using  a  laser-light- 
scattering  system  failed  to  define  the  effect  of  that  parameter  in  relation 
to  normal  experimental  error  as  a  result  of  problems  with  the  test  system 
itself  [20]. 

To  compare  the  results  of  Figure  8  with  the  available  noise  theories. 

Equation  (4)  was  used  to  describe  R^.  Here  sy  is  assumed  to  be  constant, 

which  seems  reasonable  after  inspecting  the  pressure  distribution  curve 

of  Figure  2  for  -o/C  between  0.8  and  D.9.  In  Equation  (4),  C  was 

pmin  p 

set  equal  to  C  =  -1.0,  while  the  high-speed  films  and  viedotapes 
pmin 

showed  0^  =  1/3  to  adequately  represent  the  average  pressure  experienced 

by  the  bubbles  during  collapse.  Employing  these  values  in  the  classical 

incompressible  noise  theory  given  by  Equation  (8),  a  comparison  is  made 

-8  4 

in  Figure  9  with  the  empirical  o  ’  dependence.  For  the  range  of  a 

observed  in  this  study,  the  combination  of  Equations  (4)  and  (8)  appear 

to  provide  a  reasonable  model  for  the  functional  dependence  on  a.  However, 

2  4 

Equation  (8)  also  predicts  the  noise  per  collapse  to  vary  as  Vw  ’  for 

constant  -a/C  ,  a  trend  which  was  not  revealed  in  Figure  8. 

Pmin 

Assuming  a  maximum  value  for  P  of  0.001  atm  and  using  estimates 

8 

of  R^  from  the  video  tapes,  both  Equations  (11)  and  (12),  after  insertion 
in  Equation  (15),  show  cutoff  frequencies  lying  above  100  kHz  in  the 
shock  wave  spectrum.  Hence,  Equation  (16)  should  be  a  valid  approximation 
for  the  two  shock  wave  noise  theories.  Setting  =  -1.0  again  in  Equation  (4) 
and  inserting  it  into  Equation  (17),  the  two  theories  are  compared  with 
the  measured  trend  in  Figure  9.  These  models  also  show  a  reasonable 
correlation  to  the  experimental  trends  of  the  noise  with  o.  Furthermore, 
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the  results  indicating  that  L  is  independent  of  for  constant  -a/C 

^min 

is  correctly  implied  by  using  the  combination  of  Equations  (4)  and  (17). 

Quite  possibly,  Equation  (4)  is  inadequate  for  modelling  the  actual 

maximum  bubble  radii  observed  in  this  investigation.  Using  the  videotapes, 

attempts  were  made  to  investigate  the  dependence  of  on  a  and  V^. 

Unfortunately,  valid  results  were  hampered  by  difficulties  in  obtaining 

radius  measurements  for  any  significant  range  of  o.  These  difficulties 

were  due  to  the  limited  resolution  capabilities  of  the  video  system  in 

addition  to  the  natural  distortion  of  the  bubbles  in  the  flow  field. 

Because  is  a  critical  parameter  in  each  of  the  noise  theories,  it 

is  thus  difficult  to  arrive  at  a  valid  appraisal  of  any  of  the  models 

given  the  uncertainty  of  this  parameter  as  measured  in  this  investigation. 

Since  no  other  known  studies  provide  this  necessary  data.  Equation  (4) 

stands  as  the  only  functional  approximation  of  1^,  despite  its  inadequacies. 

Energy  density  spectra  of  the  cavitation  noise  were  found  to  be  flat 

from  10  to  100  kHz  for  all  values  of  a  and  V  .  A  typical  example,  after 

corrections  were  made  for  the  hydrophone  response  and  tunnel  resonances, 

is  shown  in  Figure  10.  The  spectra  will  increase  uniformly  in 

level  for  decreasing  o.  Morozov  [18]  proved  that  the  shape  of  the  cavitation 

noise  spectrum  is  identical  to  that  of  a  single  bubble  collapse  when  the 

occurrence  of  noise  bursts  follows  a  Poisson  distribution.  Studies  of 

the  tape  recordings  indeed  showed  such  a  distribution  to  exist. 

Thus,  the  low-frequency  roll-off  predicted  by  Equation  (8)  was  not 

observed  in  the  spectra.  With  bubble  life  times  typically  around  4  msec, 

a  spectral  peak  might  be  expected  near  .25  kHz  such  that  a  roll-off  would 

then  be  anticipated  above  10  kHz.  Mellen  [11],  and  more  recently  Blake 

_2 

et.  al.  [2],  reported  an  f  trend  at  frequencies  immediately  above  the 
spectral  peak  which  was  located  near  the  reciprocal  of  the  bubble  life 
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time.  However,  the  spectra  obtained  by  Blake  et.  al.  appear  to  rise 
again  to  a  new  plateau  following  the  initial  roll-off.  Moreover,  both 
Barker  [1]  and  Mellen  [11]  reported  spectra  which  levelled  off  after 
the  low-frequency  roll-off  to  a  plateau  extending  from  10  to  100  kHz. 

Barker  also  noted  that  the  noise  level  between  10  and  100  kHz  was  essentially 
independent  of  velocity  when  a  was  held  constant,  even  though  it  increased 
as  V  ^  between  1  and  10  kHz.  Although  the  rate  of  collapse  of  individual 

oo 

bubbles  was  not  measured  in  that  study,  the  data  were  obtained  for  incipient 
conditions  where  the  collapse  rates  should  have  been  comparable.  Neverthe¬ 
less,  Barker's  results  are  also  supported  by  Figure  7. 

The  presence  of  a  high-frequency  spectral  plateau  implies  the  existence 
of  pressure  pulses,  the  time  constants  of  which  are  short  compared  to 
the  reciprocal  of  the  upper  cutoff  frequency  with  spectra  described  by 
formulas  such  as  Equation  (14).  That  the  noise  does  not  increase  with 
V  for  constant  -a/C  in  the  high-frequency  regime  seems  to  support 

oo  p 

min 

the  analysis  of  Baiter  and  Esipov  and  Naugol'nykh  if  one  assumes  that 
the  maximum  bubble  radius  depends  only  on  a.  Also,  the  a  dependence 
of  the  noise  per  collapse  is  crudely  predicted  by  the  shock  wave  models. 

On  the  other  hand.  Barker's  results  suggest  that  different  mechanisms 
generate  noise  at  low  and  high  frequencies.  Both  Mellen  and  Blake  et. 
al.  observed  spectral  behavior  supporting  the  classical  model  at  low 
frequencies.  In  addition,  the  increased  velocity  dependence  at  low  frequencies 
reported  by  Barker  is  predicted  by  a  comparison  of  Equation  (8)  to  Equation 
(17). 

CLOSING  REMARKS 

An  attempt  was  made  to  suggest  appropriate  scaling  laws  for  the 
high-frequency  region  of  travelling-bubble  cavitation  noise.  Although 
the  literature  abounds  with  numerical  schemes  yielding  the  sound 
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radiated  from  a  collapsing  bubble,  they  are  not  well  suited  for  practical 
application  to  flow-induced  cavitation.  The  pressure  field  described  in 
this  study  is  typical  of  those  associated  with  many  headforms  in  that  the 
cavitation  nuclei  experience  an  abrupt  low-pressure  impulse  that  initiates 
the  life  cycle  of  the  bubble.  Thus,  these  results  should  prove  helpful  in 
predicting  the  noise  generated  by  similar  headforms  used  in  the  travelling- 
bubble  regime.  Comparisons  with  the  available  noise  models  would  be 
improved  with  a  study  of  how  the  maximum  bubble  radii  depend  on  the  flow 
parameters  such  as  P  .  Also,  phenomena  involved  in  the  final  stage  of 
collapse  require  better  understanding,  as  this  is  the  point  where  sound 
emission  occurs.  For  instance,  thermodynamic  effects  should  be  investigated, 
because  such  an  omission  makes  even  the  models  employing  numerical  solutions 
incomplete. 
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